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For Rayleigh-Benard convection in a square cell with a fluid of Prandtl number one, we report 
experimental results on the transition between a stationary pattern of ideal straight rolls (ISR) and 
the spatiotemporal chaotic state of spiral defect chaos (SDC). In contrast to experiments in circular 
geometries, we found ISR states below a particular value of the control parameter and SDC states 
above this value. By characterizing the pattern with a global measure, the pattern entropy, we 
found that the transition from SDC to ISR showed similarities to phase transitions in equilibrium 
finite-size systems. 
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In recent years scientists in a number of fields have 
shown great interest in dissipative pattern-forming sys- 
tems of large spatial extent. These systems often dis- 
play particularly intriguing states of persistent time- 
dependent behavior termed spatiotemporal chaos [^|-Q. 
One important question regarding these states is whether 
ideas from statistical mechanics may be useful for un- 
derstanding the complex spatiotemporal behavior. In 
this Letter we present experimental data showing that a 
transition between a spatiotemporal chaotic state and a 
well-ordered time-independent state of Rayleigh-Benard 
convection is similar to phase transitions in finite-size 
equilibrium systems [Q. 

In a Rayleigh-Benard convection experiment, a hori- 
zontal fiuid layer of height d is confined between a top 
plate of uniform temperature Ttop and a bottom plate of 
uniform temperature Tbottom = ^top + ^T. For AT less 
than a critical value ATc, the fiuid is stationary; however, 
when AT > ATc, the stationary state is unstable and a 
pattern of convection rolls with wavenumber k ~ ir/d de- 
velops [y] . For systems of infinite extent in the horizontal 
directions it was shown theoretically that an ideal pat- 
tern of straight rolls (ISR) of wavenumber k is stable in 
a restricted region of the e — k parameter space g , where 
6 = {AT/ ATc — 1) is the reduced control parameter. Ex- 
periments with argon gas in a rectangular cell of small 
aspect ratio F = (length/2) /d ^ 15 [0] showed reasonable 
agreement with these predictions. However, for a larger 
circular cell of radial aspect ratio F = r/d = 76, Morris 
et al. found a state of spatiotemporal chaos within the 
boundaries of the stable region of the e — k parameter 
space [g[. This state, termed spiral defect chaos (SDC) 
and shown in Fig. [^(a), is characterized by a complex 
dynamics involving targets, rotating spirals, dislocations, 
disclinations, and grain boundaries. However, as shown 
in Fig. 0(a), SDC also shows narrow regions of slightly 
curved ISR. In this Letter we argue that it is the growth 
of these regions that characterizes the transition from 




b 



FIG. 1. (a) Spiral defect chaos and (b) ideal rolls at 
e = 0.92. Warm upflow corresponds to dark and cold down- 
flow to bright. 



SDC to ISR. 

Since the discovery of SDC an unresolved question was 
whether the ISR predicted to be stable for an infinite sys- 
tem could be realized in a large aspect ratio experiment. 
Recently, Cakmur et al. used specially prepared systems 
and showed that ISR were indeed stable where predicted, 
and that both SDC states and ISR states existed for the 
same value of e P] . Fig. [l| shows two shadowgraph images 
illustrating this bistability. 

Experiments in circular cells with fiuids of small 
Prandtl number have suggested that above a particular 
value of e, SDC is the state chosen for almost all initial 
conditions |]1Q|-|12|. Just below this value of e a time- 
dependent, apparently chaotic state without spirals was 
found. Earlier experimental and theoretical investiga- 
tions of the latter state |l[ led to the conclusion that the 
tendency of rolls to align perpendicularly to the bound- 
aries necessarily results in roll curvature which, due to 
large-scale fiows, leads to the observed persistent pattern 
dynamics. Motivated by these earlier observations, inves- 
tigations of SDC were mostly limited to circular contain- 
ers. Almost no experiments were conducted in rectangu- 
lar cells [|l2| , although pioneering experiments by Gollub 
and coworkers [13] had shown that ISR appeared after a 





FIG. 2. Shadowgraph images of the convective pattern for 
a control parameter jump from below onset (e < 0) to above 
onset (e = 0.233) at times (a) - 45rT, (b) 229tt, (c) 458tt, 
and (d) 22689tt = 9.1rh, the stationary pattern. 



long transient in a small aspect ratio rectangular cell. In 
the remainder of this Letter, we present quantitative ex- 
perimental data from a large aspect ratio convection cell 
of square geometry showing that the transition between 
SDC and ISR exhibits similarities to equilibrium phase 
transitions in finite-size systems. 

The experiment was conducted in a square cell with 
aspect ratio F = {L/2)/d = 50, where L is the length on 
a side. The fluid was CO2 gas at a pressure of (41.593 ± 
0.007) bar and Prandtl number a = tt/tj^ = 1.1, where 
Ti, is the vertical viscous timescale, tt = (P/k, = 2.6s 
is the vertical thermal diffusion timescale, and Hi is the 
thermal diffusivity [Q . The experimental setup was sim- 
ilar to that described in Ref. ^^. The cell's circular top 
and bottom sapphire plates were 1 cm thick, were spaced 
(623 ± 4) /im apart, and were parallel to within ±3/im 
over the 10 cm diameter. The top plate temperature was 
set at (24.00±0.05)°C by a circulating water bath and the 
bottom plate was heated by an electric film heater. Both 
temperatures were regulated to ±0.3 mK, and the pres- 
sure was regulated to ±5 x 10~^ bar. For this experimen- 
tal situation the parameters were sufficiently tempera- 
ture independent so the Boussinesq approximation could 
be applied p4|Jl5| ] . The bottom sapphire plate was coated 
with aluminum to allow the visualization of the pattern 
from above using the shadowgraph technique Q . Eight 
circular paper sheets were placed between the top and 
bottom plates, and a square of size lOOd was cut out of 
the center of the circular sheets to provide the boundary 
of the convection cell [16|. The measured onset of con- 
vection was (2.03 ± 0.02) K, which is in good agreement 
with the theoretical prediction of (2.04 ± 0.02) K. 



FIG. 3. Time evolution of the pattern for e — 0.58 at 
times (a) 13.70r/^, (b) 13.80t/,, (c) 14.07r^, and (d) 28.96r^, 
the stationary pattern. 



For each value of the control parameter studied, the 
system was initialized by a jump from below onset (e < 0) 
to above onset (e > 0). After this jump, straight rolls ini- 
tially formed near the sidewalls while a random pattern 
appeared in the middle of the cell. For e ^ 0.53, the 
pattern coarsened over time and developed after a tran- 
sient into ISR with a few stray defects. An example of 
the pattern evolution is shown in Fig. g. For e ^ 0.58 
the pattern did not order into ISR during observation 
times of at least 20t/i, where th = V^tr is the horizon- 
tal thermal diffusion time. In the intermediate regime, 
0.53 ^ e ^ 0.58, the behavior was more complicated with 
large patches of almost ISR growing and competing with 
large patches of SDC. Sometimes, often after long tran- 
sients, a patch of ISR grew and filled the entire system. 
As can be seen in Fig. ||, once a patch of ISR connected 
two opposing sidewalls it grew via the propagation of two 
almost flat fronts. The behavior described above is rem- 
iniscent of that of an equilibrium finite-size system near 
a phase transition ^ . 

We note that the initial pattern evolution showed the 
brief appearance of spirals for values of the control pa- 
rameter as low as e ~ 0.2 (as shown in Fig. [^(b)). This 
observation suggests that the bistability of ISR and SDC 
may extend to lower control parameter values than those 
reported here. This conjecture would be in agreement 
with earlier observations in a larger aspect ratio circular 
cell (F = 78) in which SDC was found above e ~ 0.22 



It would also be 
|ll|,|l2|,0 and simu- 



for similar experimental conditions 

consistent with both experiments 

lations Pq] in which it was found that the onset of SDC 

decreased with increasing system size. 

To quantify the transition from SDC to ISR, we cal- 




FIG. 4. Correlation length ^ in the spatiotemporal chaotic 
regime. 

culated several statistics in the vicinity of the transition. 
On the SDC side of the transition, we measured the cor- 
relation length ^ from the exponential decay of the auto- 
correlation function: 

C(Af) = {{u{x^Ax,t) - {u)) {u{x,t) - {u)))^^u (1) 

where u{x^ t) is the intensity of the shadowgraph picture. 
The autocorrelation function is obtained from the Fourier 
transform of the power spectrum averaged over the du- 
ration of each experimental run [Q. As shown in Fig. |^, 
the correlation length appears to diverge at a finite e. 
The solid line in Fig. | is given by ^ = 1.9d(e-0.57)-°-^2. 
This result is surprising since in an experiment with cir- 
cular geometry the correlation length was found to di- 
verge as ^ oc e~^"^^ ]q,[l4. It should be noted, however, 
that in the same circular experiment [Q the measured 
data for the correlation time was consistent with a diver- 
gence at finite e, with an apparent transition between two 
chaotic states. In our experiment, however, we observed 
a transition between a spatiotemporal chaotic state and 
a stationary state. These apparent discrepancies may be 
explained by the incompatibility of the circular geometry 
with ISR, as discussed above. 

To further quantify the transition between SDC and 
ISR, we have calculated the spectral pattern entropy 
[|2Q| , pi[ . This quantity is defined through a spectral dis- 
tribution function: 



p{k,t) = 



\^{k,t)\^ 



jdPk\^f{k,t)[^ 



(2) 



where ^(k^t) is the Fourier transform of the two dimen- 
sional pattern at time t. Using this distribution function, 
the spectral pattern entropy S{t) can be defined as: 



S{t) = - / d?kp(k,t)\np{k,t). 



(3) 



S{t) measures the disorder in the pattern. For example, 
if the pattern is ideal (only one mode is excited) 5 = 
and otherwise 5 > 0. 

In Fig. pi the time evolution of the pattern entropy is 
shown for e = 0.554. The patterns at the marked places 
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FIG. 5. Time evolution of the pattern entropy S{t) and 
the pattern for e = 0.554 at times (a) 14.52r/i, (b) 20.27rh, 
(c) 30.95T/,, and (d) 32.32r/,. 



show a clear correlation with the value of the pattern en- 
tropy. As can be seen in Figs. |^(a) and ^(c), patterns 
with disordered or curved regions lead to large values of 
S(t)^ while patterns with straight regions, such as those 
in Figs, ^(b) and ||(d), lead to small values of S(t). The 
latter two patterns also show that during the evolution, 
well-ordered regions were typically aligned either diag- 
onally or perpendicularly to one of the sidewalls of the 
cell; it appeared that the pattern was probing the sys- 
tem's symmetries. 

The experimental data are summarized in Fig. g [ p2[ . 
Fig. ||(a) shows the temporal average of the pattern en- 
tropy {S)t as a function of e. As the transition to ISR 
is approached from above, {S)t shows a sharp decrease. 
Fig. ^(b) shows the variance cf{S) over the same range of 
e. As the transition is approached, (t{S) displays a sharp 
increase. The location of the sharp changes is consistent 
with the divergence of the correlation length in Fig. 0. 
For larger e the variance cf{S) approaches a small value, 
suggesting that the system consists of many independent, 
fluctuating subsystems. Again, the behavior shown in 
Fig. ^ appears to be similar to that of equilibrium phase 
transitions Q]. 

For a fluid of Prandtl number a = 1.1 in a square con- 
vection cell of large aspect ratio, we have found a compe- 
tition between spiral defect chaos and ideal straight roll 




FIG. 6. (a) The time- averaged pattern entropy {S)t and 
(b) variance cr(S) as a function of e. The inset in (a) shows 
data for both SDC and ISR. 



states. This competition is particularly evident in the 
interplay of disordered regions and straight roll regions 
during the evolution of the SDC states. Using spatial 
correlation lengths and statistics based on the spectral 
pattern entropy, we have provided quantitative evidence 
for a transition between SDC and ISR at a finite value of 
e. This transition shows many similarities to transitions 
in finite-size equilibrium systems, including an intrigu- 
ing intermittent behavior near the transition. We are 
currently investigating whether experiments in cells of 
different aspect ratios can be used to predict the onset of 
SDC in infinitely extended systems. 
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